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• Motivation and background 

• Description of peak-seeking algorithm 

• Implementation on F/A-18 

• Performance data flight 

• Simulation results 
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Introduction 



• US domestic flights in 2011: 

- 12.1 billion gallons of fuel 

- 114.6 million metric tons of C0 2 equivalent 

• NASA's Environmentally Responsible Aviation project 

- Mitigate the impact of aviation on environment 

- Reduce fuel consumption, emissions, and noise 

• Concept presented here: 

- Reduce drag in cruise by altering the trim configuration, 
applicable to many types of aircraft 
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Background 



• Existing Trim Methods 

- Often scheduled with flight condition 

- Based on a priori information (analytic, wind tunnel, flight data) 

- Differences between models and reality may degrade performance 

• Off nominal flight conditions, lifetime variations, manufacturing 
differences, external modifications or stores, etc... 

• Real-time optimization methods 

- Adaptive Performance Optimization 

• Drag reduction on L-1011 by use of symmetric aileron, (Gilyard et al.) 

- Formation flight 

• Position optimization (Ryan and Speyer) 

• Spanwise lift distribution optimization (Hanson and Ryan) 

- Trim optimization 

• Drag reduction by use of single trailing edge surface group on X-48, in 
simulation (Griffin et al) 
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Approach 


• Real-time optimization of trim configuration to reduce 
drag 

• Use any number of control effectors 

• Utilize onboard measurements of performance, which 
may be noisy 
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Performance Measurement, / 


Peak-seeking Scheme (simplified for 1 effector) 
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Peak-seeking algorithm 
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Technical Formulation: Performance Function 


Performance Function (Taylor series): 

/(%) ~ f(Xk- 1) + £>l (x k - - x fe _ a ) 


Assuming the performance function can be treated as linear at any control surface 
position and expanding to include any number of control effectors, n, gives: 
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F and x are measureable, b k is unknown and to be estimated, and since F and x are 
noisy and F varies with x, a time-varying Kalman Filter is an appropriate choice for 
an estimator. The states of the Kalman filter are define as the gradient vector: 
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Technical Formulation: Kalman Filter 


Measurement equations are expanded to include multiple previous measurements, M: 



, -se 
1 

*1—1 

1 

T 

X 2fc_i *2^ ■■■ X nfc-1 

t> 

**1 

II 

/fe_ 2 ) - /Or*) 

II 

a? 

X lk- 2 ~ X lk X 2k-2 ~ ■■■ X n k -2 ~ X nk 

■ ■ -i i 

: : ** : 


- f(x k - M ) - f(x k ). 


X lk-M - X lk X 2 k-M ~ X 2 k ■" X n k -M - X nk- 


Kalman filter measurement equation: Kalman filter process equation: 

AF^ = + Ujz ~ Cfe-i + w k 

where v k , w k are Gaussian white-noise with covariance matrices R k and Q k respectively 


A standard linear time varying Kalman filter is 
then implemented as follows: 
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Persistent Excitation and Initial Excitation 


• Persistent Excitation 

- Addition to commanded 
surface positions that is 
helical about the trajectory 



• Initial Excitation 

- M points around a 
circle/sphere centered at 
the initial condition 
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F/A-18 : NASA 853 


• Modified F/A-18 Aircraft - Research flight control computers 

• Nonlinear Dynamic Inversion inner loop control laws 

• Autopilots: 


• Altitude Hold 

• Airspeed Hold 

• Wing Leveler 
algorithm adds biases 

• Symmetric aileron 

• Trailing-edge flaps 

• Leading-edge flaps 




Inboard/outboard leading 
edge flaps (ganged together) 


Ailerons 
Trailing edge flaps 


Stabilators 
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Performance Data Flight 


• Early in development an opportunity was 
presented to collect performance data during 
another research activity's flight. 

• Commanded 80 test points with combinations of 
leading edge flaps, trailing edge flaps, and 
symmetric ailerons and recorded resulting fuel 
flow over >30sec per pt. 


• Evaluated at a single flight condition of 25,000ft, 
240 KCAS 
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Performance Model 


• Developed a new plant model for simulation testing. 

• Polynomial fit to flight data across 3 axes 
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Fuel flow, % change from baseline 



Performance Model 


• More detailed data set collected for trailing edge 
flaps vs symmetric ailerons, leading edge at 5 deg 

— Spanwise lift distribution control 






Baseline 

- Trailing edge flaps, 5 to 6 deg 

- Symmetric ailerons, 0 deg 

Minimum, -2.3% 

- Trailing edge flaps, 3 deg 

- Symmetric ailerons, 5 deg 
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Fuel flow- % change from baseline 


Fuel flow 


Noise Model 


• Generated a noise model for simulation, added 
onto output from new performance plant model 
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New plant model for simulation 



• Using new plant model in simulation, peak 
seeking controller was evaluated and tuned 

• Tuning variables: 

- Gain applied to gradient, "controller gain" 

- M, number of previous measurements used by 
Kalman Filter 

- R and fuel flow filter time constant, tuned for 
signal noise 

- Q, Kalman filter process covariance 
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Trailing-edgs flaps, deg 


Gain Tuning 
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M, previous measurements 
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Q, R, Fuel flow filter 



• Filter on fuel flow time constant and R matrix 

- filter to reduce noise on signal going into Kalman 
filter, adjust R accordingly 

• Q matrix, process covariance, tuned through 
Monte Carlo type simulation 
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Final Tuned Parameters 


Parameter 

Value 

Gain 

-105 

M 

5 for 2 effectors 
7 for 3 effectors 

Fuel flow filter 
time average 

20 s 

R 

1.85 2 I 

Q 

1.98 2 I 
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Simulation Results - 2 effector 


Starting from 4 different positions, algorithm 
converges around -2% 
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Symmetric ailerons, deg 
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Fuel flow, %change from baseline 


Simulation Results - 2 effector, case B 
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T railing-edge flap, deg Leading-edge flap, deg 


Simulation Results - 3 effector 


3 effector test, converges to -2.5% 
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August 13, 2013 symmetric aileron, deg 


- Unfiltered Fuel Flow 
Bias Inserted 

- Algorithm Engagement 



Conclusions 


• Peak-seeking algorithm has potential to reduce 
fuel consumption on wide variety of aircraft types 


• Can easily be implemented into existing control 
structure (assuming ability to actuate multiple 
effectors, and digital control) 


• Algorithm was subsequently flown on 5 flights 
accumulating about 5 hours worth of test data 

- Results will be presented tomorrow at 5:30pm (Salon J) 
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